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ABSTRACT 

We study the nonperturbative functions in the Colhns-Soper-Sterman resum- 
mation formahsm by examing Drell-Yan data in both fixed target and coUider 
experiments and then predict the transverse momentum distributions of the 
and bosons at the Tevatron. Our results differ from that in the hterature and 
agree better with pubhshed CDF data. Using statistical arguments, we find a 1 fb~^ 
luminosity at the Fermilab Tevatron should be able to provide useful constraints 
on the nonperturbative functions. 



1. Introduction 



The measurement of the mass (M^) of the W boson is very important in 
testing the Standard Model (SM), which has proven to be extremely successful. 
This measurement is currently done at the Tevatron by both the CDF and the DO 
groups. In the framework of the SM, the mass of the Higgs boson can be known 
to within a couple of hundred GeV if M\y is measured to within 100 MeV and the 
mass of the top quark within 10 GeV. 

Since the boson decays into a charged lepton and a neutrino, it is important 
to know the kinematics of boson production for measuring Mw, in particular, 
the transverse momentum distribution of the bosons. The measurement of the 
transverse momentum of the gauge bosons V (either W^, or virtual photon in 
the Drell-Yan process) also provides a test of QCD in processes involving two-scale 
problems which are beyond the usual framework of perturbative calculations in the 
expansion of the strong coupling ag- The measurement of the rapidity asymmetry 
in the production of bosons provides a handle on the ujd ratio in parton 
distribution functions such that quantities like V-\^ /Vz can be determined to good 
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precision. 

To obtain the kinematics of the gauge boson we must resum the multiple soft 
and colhnear gluon effects. We adopt the CoUins-Soper-Sterman (CSS) formal- 
ism'^'*' to resum these multiple gluon effects and closely follow the notation used in 
Ref. 4. In this paper we study the nonperturbative part of this resummation for- 
malism and show the importance of its contribution to the transverse momentum 
distribution of the gauge bosons. 

To describe the kinematics of the gauge boson for all transverse momenta Qt-, 
we use the rcsummed formula differential in the mass [yj Q"^), rapidity (y), and Qt 
of the boson Vj"'' 



X exp 



In 



Q 



^ J hqib) - hj/A(xA, b) - hyB(xB, b) 



+ Y{Qt;Q,xa,xb) >, 



(1.1) 

where ctq provides the process-dependent normalization and W and Y are derived 



from pertrubative calculations. 



[5] 



2 



The functions /ij/^ and which carry a flavor dependence as well as their 

respective dependence on the momentum fractions xa and xb, handle the nonper- 
turbative behaviour at large b along with hq. These functions are to be obtained 
by a fit to data, subject to the constraint that they must vanish when 6 — > 0. The 
constant Qq, in Eq. (1.1), is completely arbitrary. 

Among the nonperturbative functions used previously in the CSS formalism 
are those of Davies, Stirling and Webber !^'^' They selected the functional form that 
provided a gaussian smearing of the transverse momentum, 

hqib) = and hj/A{xA: b) + hj/s{xB: b) = ^16^, (1.2) 

where gi,g2 are phenomenological constants. For simplicity the nonperturbative 
functions are assumed to be independent of flavor and momentum fractions xa, 
Xb- Using the parton distribution functions of Duke and Owens, their favored 
values after fitting to and R209'"" Drell-Yan data with Qq = 2 GeV and 

bmax — 0.5 GeV~^ were 

gi = 0.15 GeV^ and §2 = 0.4 GeV^. (1.3) 

The interest at that time was the production of bosons in hadron coUisons at 
^/s = 540 GeV. 

2. A New Study 

For this study we chose Qq — 1.6 GeV and bmax = 

0.5 GeV~^ with a different 

functional form for the nonperturbative part, 

hqib) = g2b^ and hj/A{xA, b) + hj/B{xB, b) = gib[b + g3ln{100xAXB)]- 

(2.1) 

where gi,g2,g3 are phenomenological constants. Since we want to investigate the 
production of gauge bosons at the Fermilab Tevatron with ^/s — 1.8 TeV, our 
kinematic region includes values of r = xax b that are significantly lower than those 
relevant to Ref. 6. For this reason we refit for the nonperturbative functions using 
lower mass ranges in the R209 {pp ^~^fi^+X at ^/s = 62 GcV) and E288 {pN 
fi^fi~ + X at ^/s = 27 A GeV) data in addition to the data from CDF {pp collisions 
at ^/s = 1.8 TeV) for Z'^ production. These lower r values used from the R209 and 
E288 data, however, are still roughly a factor of two larger than the typical r value 
of and Z° physics at the Tevatron. To improve these nonperturbative functions 
in the CSS formalism, we postulate the xa,xb dependence for the nonperturbative 
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functions hjj^ and h-jj^ as in Eq. (2.1). Our choice of \yi{xaxb) is inspired by the 
fact that the average transverse momentum of the Drell-Yan pair grows slowly with 

mi 

T. When better statistics becomes available, these nonperturbative functional 
forms also can be extracted for the lower r kinematics through the study of Drell- 
Yan pairs and boson production at the Tevatron collider. 

For boson studies at the Tevatron, it would be excellent if one had very 
high statistics for the data, since this would enable a determination of the 
nonperturbative functions in a kinematic region that overlaps greatly with the 
kinematics for boson production, but such high statistics data do not exist. 
Using the R209 data over the mass range 5 < Q < 8GeV in conjucntion with 
the Z^ boson data from CDF, g2 was determined. We note that g2 is associated 
with the ln(2^) factor in the nonperturbative functions. With that g2, the same 
R209 data were taken in combination with the E288 data over 6 < Q < 8 GcV 
to obtain an xa^xb dependence by determining gi and g^. Using the CTEQ2M 
parton distribution functions (PDFs) with Qq = 1.6 GeV and bmax = 

0.5GeV-\ 

the nonperturbative parameters are 

gi = 0.1lig04 GeV2, g^ = 0.58^^^ GeV^, ^3 = -l-Stg } GeV-^. 

(2.2) 

We now proceed to discuss the comparison of the Qt distributions obtained 
from the previous values of gi,g2 with that given by gi,g2,g3- For the high r 
kinematics of the R209 data with 11 < Q < 25 GeV, the da/dQ^ distribution 
given by the nonperturbative form of Eq. (1.3), where we naivley carry over the g 
to the CTEQ2M PDF, lies very close to the results provided by using the CTEQ2M 
PDF with Eq. (2.2). The results provided by Eq. (2.2) also are consistent with 
the E288 Drell-Yan data, whose typical r values range from about | to |. For 
the 5 < Q < 8 GeV range in the R209 data, where the r values probed are much 
smaller and more relevant to W"^ and Z'^ boson physics at the Fermilab Tevatron, 
Fig. 1 shows that the old form of Eq. (1.3) with the CTEQ2M PDF is quite different 
from the result of Eq. (2.2), deviating from the experimental result. This is a clear 
demonstration that for precision measurements (like the determination of M]y or 
r^) or simply for theoretical consistency, it is necessary, particularly when entering 
a new kinematic region or changing the PDF used, to have a nonperturbative 
paramctrization consistent with that PDF and kinematics. In Fig. 1 we also show 
the result obtained when the favored values of Ref. 6 are used in conjuction with 
the HMRSB PDF, ' ' which has been used in previous publications. This too does 
not agree with the data, clearly discriminating between the two nonperturbative 
functions of Eq. (1.3) and Eq. (2.2). The former parametrization is invahd in this 
kinematic region independent of the PDF used. 
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Having a better nonperturbative dependence for describing the low Qt kine- 
matics of the W^, boson production at the Fermilab Tevatron, we can see if 
there is any change from previous expectations. In Fig. 2a (2b) we display three 
calculations against the results for (W^^) production obtained by the CDF col- 
laboration. (We use Mw = 80 GeV and Mz = 91.17 GeV.) The dashed (dashed- 
dotted) curve is where we took the previous values of Eq. (1.3) with the CTEQ2M 
(HMRSB) PDF, while the solid curve represents the results of the new fit with the 
CTEQ2M PDF. Simply comparing the new set with the old set as applied to the 
CTEQ2M PDF, it is apparent that the peak of the Qt spectrum has dropped and 
shifted to a higher value. The HMRSB result has a lower peak height, but 
the peak remains at the lower Qt, just as the dashed curve. It is not due to the 
comparison of these theoretical results with the low statistics data that favor one 
result over the other, rather, it is the inability of the Eq. (1.3) to account for the 
Drell-Yan results with the more pertinent r kinematics that decides. Recall that 
at very low Qt the da/dQT experimental results'"' for the boson were below 
theoretical expectations with the old fits and that as the Qt rose the data quickly 
gave experimental results above theoretical expectations!^^' The new fit not only 
shifts the peak in the da/dQT distribution to higher Qt, it also makes the peak 
broader, thereby yielding an improvement between the data and the theory in the 
region of 5 < Qt < 20 GeV. In Fig. 3a (3b) we show the Qt distribution at y = 
for the Z° (W^) boson relevant to the CERN collider!'" 

We observe that the results with the old parameters are similar between the 
CTEQ2M and HMRSB PDFs for the r region of Fig. 1, yet as we enter the lower r 
kinematics of the CDF data, the difference appears to grow. In this regard, we note 
that the CTEQ2M PDF fits with the low-x F2 data from HERA, while HMRSB 
does not. Using the new nonperturbative parameters with the MRSD— ' PDFs''*' 
gave results that differed only negligibly when compared against the CTEQ2M 
results for Qt < 20 GeV. 
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3. Measuring the Nonperturbative 
functions at the Fermilab Tevatron 



Measuring the nonperturbative functions in the CSS resummation formahsm 
provides information about the nonperturbative nature of QCD theory. For in- 
stance, hq^b) in Eq. (1.2) might have an operator definition as the vacuum expec- 
tation of the gluon condensate!^^' 

We explore the necessary luminosity and energy of the upgraded Tevatron to 
perform this kind of study. To isolate the problem interested, we assume that 
the parton distribution functions are known in the relevant kinematic region for 
producing and bosons at the Tevatron. In Fig. 4 we show the transverse 
momentum distributions for the and bosons produced at a pp collider with 

= 1.8, 3.5 TeV using the CTEQ2M PDFs. 

For this study we chose Qq = 1.6 GeV, fixed = —1.5 GeV~^ and pretended 
that nature wants (71 = 0.11 GeV^ and g2 = 0.58 GeV^. To estimate the accuracy 
for the measurement of the nonperturbative parameters as a function of luminosity, 
fake data was generated for the da/dpT distribution describing the production of 
Z^ bosons, where a branching fraction of 0.06 was included to focus only on the 
Z — > e~e^ and Z IJT ji^ decay channels. This fake data was then fit to theory 
using MINUIT'^°'to both determine the best values for the parameters gi and g2 
and to estimate the errors. 

To create the fake data for the da/dpT distribution, a bin width of ApT = 
1 GeV was assumed. Given the theoretical value for da/dpj^ and the luminosity 
{£), a statistical error was evaluated for each value of da/dpT, 



Assigning a detector uncertainty equivalent to the statistical error for the purposes 
of this estimate, the error for each of the fake data points in da/dpT becomes 



Taking eight values at pt = 1, 2, . . . , 8 GeV, each point of the da/dpT theory 
was randomized according to a gaussian distribution about its theoretical value 
using a width of e. This provided a fake data sample to use for the fit of the 
nonperturbative functions. 




(3.1) 




(3.2) 
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When the fake data was fit to the theoretical values, the errors provided by 
MINUIT demonstrated that gi can be known to within about 6% (20%, 50%) and 
g2 can be known to within about 2% (4%, 15%) at the 95% confidence level given 
a luminosity of 10 fb"^ (1 fb"\0.1 fb"^). 

The accuracy for the measurement of the nonperturbative functions in the 
CSS resummation formalism is relatively indifferent to the energy upgrade of the 
Tevatron (within a factor of two for 1.8 TeV and 3.5 TeV energies) when compared 
against an orders of magnitude increase in the luminosity, indicating that the 
luminosity of the machine is crucial. In this analysis, we only assume Z e^e~ 
and Z — > data. If the mass of the W boson is known, one can include the 

data from + neutrino (for Z = e or //) to further improve the measurement 

by about a factor of two. 

Although this study is merely a theoretical exercise, the nonperturbative func- 
tional forms as described in Eq. (2.1) might need to be revised when better data 
become available. It eventually may prove that the flavor dependence cannot be 
ignored. Nevertheless, our point is that data from Drell-Yan experiments and the 
production of and Z° bosons at hadron coUiders and fixed target experiments 
like those at Fermilab will be instrumental in providing bounds on the nonpertur- 
bative structure at low Qt- 

4. Conclusion 

For precision measurements (like the determination of M^y, F^i/, or the charge 
asymmetries for bosons) or simply for theoretical consistency, it is necessary, 
particularly when entering a new kinematic region or changing the PDF used, to 
have a nonperturbative parametrization consistent with that PDF and kinemat- 
ics. In this study we have demonstrated that in the CSS resummation formalism 
represented by Eq.(l.l), the contributions of the nonperturbative functions to the 
Qt distribution of Drell-Yan pairs and the production of and bosons at 
fixed target and hadron colliders are important. When the physics of interest {e.g., 
low Qt boson production) probes different kinematics, such as lower momentum 
fractions xa and xb-, these nonperturbative functions have to be reevaluated, just 
as the PDFs have to be updated when new data probing smaller x regions become 
available. Consistency between the data from Drell-Yan processes and the produc- 
tion of and Z^ bosons will not only support QCD theory, but also will provide 
a tool which can facilitate our pursuits in physics beyond QCD through a bet- 
ter understanding of signal and background processes in this important kinematic 
region where event rates are large. As an example, we have shown that the theoret- 
ical Qt distributions for and Z^ bosons at the Fermilab Tevatron agree better 
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with experimental data after a new fit for the nonperturbative dependence has 
been performed using Drell-Yan data proximate to the relevant kinematic regions 
in T = xaxb for and Z° boson production at the Fermilab Tevatron. 

For the precision measurement of Mw and testing QCD theory in processes in- 
volving two-scale physics, such as the W^, and Drell-Yan pair production, it is 

important to know the theoretical errors due to the nonperturbative parametriza- 
tion, the factorization scale dependence, and the parton distribution functions. All 
of these considerations are under study!^^' 



8 



Acknowledgement s 

We thank our colleagues P. Agrawal, J. Botts, G. Brandenburg, R. Brock, 
J. Collins, K. Einsweiler, J. Huston, H. Lai, J. Ng, J. Paradise, B. Pope, J. Qiu, 
D. Soper, G. Sterman, W.-K. Tung, H. Wccrts, and B. Winer. This work was 
performed in part under TNRLC grant RGFY9240 and NSF grant PHY93-09902. 



9 



REFERENCES 



1. M. Swartz, XVI International Symposium on Lepton-Photon Interactions, 
Cornell University, Ithaca, New York, 10-15 August 1993. 

2. E.L. Berger, et al, Phys. Rev. D40, 83 (1989); erratum ibid, D40, 3789 
(1989). 

3. J. Collins and D. Sopcr, Nucl. Phys. B193, 381 (1981); erratum ibid, B213, 
545 (1983); ibid, B197, 446 (1982). 

4. J. Collins, D. Soper and G. Sterman, Nucl. Phys. B250, 199 (1985). 

5. Note Eq. (1.1) represents the nonperturbative factor differently than Ref. 19, 
and we have incorporated the flavor dependent couplings into Wj. 

6. C. Davies, Ph.D. Thesis, Churchill College (1984); C. Davics and W. Stirhng, 
Nucl. Phys. B244 (1984) 337; C. Davies, B. Webber and W. Stirling, Nucl. 
Phys. B256 (1985) 413. 

7. P.B. Arnold and R.P. Kauffman, Nucl. Phys. B339 (1991) 381. 

8. D.W. Duke and J.F. Owens, Phys. Rev. D30 (1984) 49. 

9. A.S. Ito, et al., Phys. Rev. D23 (1981) 604. 

10. J. Paradiso, Ph.D. Thesis, Massachusetts Institute of Technology (1981); 
D. Antreasyan, et al, Phys. Rev. Lett. 47 (1981) 12. 

11. J. Botts, J. Huston, H. Lai, J. Morfin, J. Owens, J. Qiu, W.-K. Tung, 
H. Weerts, Michigan State University preprint MSUTH-93/17. 

12. P.N. Harriman, A.D. Martin, R.G. Roberts and W.J. Stirhng, 
Phys. Rev. D42 (1990) 798. 

13. F. Abe, et al, Phys. Rev. Lett. 66 (1991) 2951; B.L. Winer, Ph.D. Thesis 
(1991), Lawrence Berkeley Laboratory preprint LBL-30221. 

14. F. Abe, et al., Phys. Rev. Lett. 67 (1991) 2937; J.S.T. Ng, Ph.D. Thesis 
(1991), Harvard University preprint HUHEPL-12. 

15. UA2 CoUaboration (J. Ahtti et al.), Z. Phys. C47 (1990) 523. 

16. HI Collaboration (I. Abt, et al.), DESY-93-117, Aug 1993; ZEUS Collabo- 
ration (M. Derrick, et al.), Phys. Lett. B316 (1993) 412 

17. R. Brock, G. Ladinsky, W.-K. Tung, C.-P. Yuan, in preparation. 

18. A.D. Martin, W.J. Stirling, and R.G. Roberts, Phys. Lett. B306 (1993) 145; 
erratum, ibid, B309 (1993) 492. 

19. G. Sterman, private communication. 



10 



MINUIT, Application Software Group, Computing and Networks Division, 
CERN, Geneva, Switzerland. 



11 



Figure Captions 



1. Comparison of R209 Drell-Yan data with calculations using Eq. (2.2) and 
the CTEQ2M PDF or Eq. (1.3) using either the CTEQ2M or HMRSB PDFs. 

2. Comparison of (a) boson or (b) W{— + W~) boson production at 
CDF with calculations using Eq. (2.2) and the CTEQ2M PDF or Eq. (1.3) 
using either the CTEQ2M or HMRSB PDFs. 

3. Calculations as in Fig. 2, except for ^ — 630 GeV. 

4. The transverse momentum distribution for the production of (a) Z and (b) W 
bosons in pp coUisions at ^/s — 1.8, 3.5 TeV. 
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